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Nanocrystalline bulk materials are desirable for many applications as they combine 
mechanical strength and specific electronic transport properties. Our bottom up ap- 
proach starts with tailored nanoparticles. Compaction and thermal treatment are 
crucial, but usually the final stage sintering is accompanied by rapid grain growth 
which spoils nanocristallinity. For electrically conducting nanoparticles, field acti- 
vated sintering techniques overcome this problem. Small grain sizes have been main- 
tained in spite of consolidation. Nevertheless, the underlying principles, which are of 
high practical importance, have not been fully elucidated yet. 

In this combined experimental and theoretical work we show, how the develop- 
ing microstructure during sintering correlates to the percolation paths of the current 
through the powder using highly doped silicon nanoparticles as a model system. It is 
possible to achieve a nanocrystalline bulk material and a homogeneous microstruc- 
ture. For this, not only the generation of current paths due to compaction, but also 
the disintegration due to Joule heating is required. The observed density fluctuations 
on the micrometer scale are attributed to the heat profile of the simulated powder 
networks. 



I. INTRODUCTION 



A growing demand for nanocrystalline materials is driven by applications of energy ma- 
terials e.g. for lithium batteries or fuel cells^, or for thermoelectricity^ which make use 
of specific mechanical, electrical, magnetic or optical properties of the nanomaterial. The 
bottom-up approach towards nanocrystalline bulk material starts from a nanoparticle pow- 
der and compacts it in a way that preserves the specific nanofeatures. Within the last years, 
this bottom-up approach has become increasingly applicable for a growing number of chem- 
ical compositions, as new reactor concepts and strategies can supply tailored nanoparticles 
like zinc oxide, silicon, alumina, or yttrium oxide in sufficient quantities. 

The way towards a dense nanocrystalline sample is tricky 4 . A characteristic feature of 
nanopowders is their low tap density combined with high stiffness, mainly originating from 
strong interparticle van-der-Waals forces. 

The energy stored within the surface of the nanomaterial is a driving force for coarsening 
as soon as temperature is applied. The consequence is that the specific nanofeatures get 
lost. Traditional sintering methods like pressureless sintering and hot pressing require the 
samples to be kept at high temperatures for a certain time resulting in strong coarsening, so 
that the effort of synthesizing nanosized starting particles is lost during the sintering process. 
Sometimes, the over-all result is even inferior to sintering micrometer-sized particles. 

As a consequence, an old method has been rediscovered, compare review^. Field activated 
sintering techniques (FAST), also called electric current sintering or spark plasma sintering, 
have demonstrated their strength in sintering nanoparticles to dense bulk under preservation 
of the nanogranular structure, while densities near the one of the respective crystalline 
bulk material can be obtainecP. Field assisted sintering is a rather quick method, which 
makes it also attractive for production processes. For both, hot pressing and FAST, the 
compaction is done by a combination of pressing and heating. The sample is located within 
a graphite crucible between graphite dies. Those create reducing conditions. During the 
whole process, a pressure is applied on the sample, usually in the order between lOMPa 
and 100 MPa, which, by itself, is too small to compact the powder. For highly insulating 
material, the basic difference between hot pressing and FAST is the heating process. Field 
assisted sintering techniques use Joule heating applied by driving a kA-current through the 
graphite crucible, and ohmic heating increases the sample temperature from outside. This 

2 



creates a radial temperature field within the sinter-bodj^l Still, the same main mechanisms 
as observed for hot pressing, lattice diffusion and grain boundary diffusion, dominate the 
sintering process^. 

If conducting bulk material is exposed to the electric current heating, it is found that the 
temperature gradients within the sample are much higher than in non-conducting material, 
shown by temperature field simulations 7 . If particles shall be compacted using field assisted 
sintering, the intrinsic mechanism of densification becomes a different one. Especially, if 
the process is conducted in a way that the sample is electrically isolated from the crucible 
and the path of the electric current is through the sample preventing a short cut via the 
crucible. Then percolation effects inside the sample have to be considered. For a given 
current through two particles in contact, most ohmic heat is released, where the resistance 
is highest, which is generally thought to be at the interface. As a result, a partial melting 
at the surface takes place as supposed by Groza et al. | 10 | n | . After the current pulse, this 
partially molten zone will have a lower resistance — with back coupling on the current path. 

The successful utilization of this process has been reported, but a theoretical explanation 
how it works is still missing. In a combined experimental and theoretical approach, we have 
investigated the field-assisted sintering of highly doped silicon nanoparticles. Nanocrystalline 
silicori^n^l and silicon based alloys^ are promising materials for thermoelectric applications. 
Additionally, they are also a model system for a conducting, covalently bound material. It 
has been shown by several groups that nanocrystalline silicon can be densified applying field 
assisted sintering techniques with very short sintering time in the order of minutes^DSl. 

II. METHODS 

II. 1. Model of field assisted sintering 

Electrical currents through a nanopowder flow along a complex, inhomogeneous and 
anisotropic network of paths, essentially determined by the pore structure. Hence Joule 
heating is non-uniform and leads to structural changes correlated with the structure itself in 
an intricate way. These correlations between structural kinetics and local structure are an 
important aspect of field assisted sintering. They influence the outcome of the compaction 
process, as we are going to show below. 
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To investigate such statistical correlation mechanisms of field assisted sintering, a phe- 
nomenological network model may be used. Densification, the current path through the 
powder, the obtained temperature fields, and the structural evolution can be calculated. In 
this approach material properties and processing enter only via a few thermodynamic and 
electrical parameters (e.g. specific heat, melting temperature, resistivity, sample cooling, 
pressure). Comparison with microstructures obtained experimentally for different sintering 
parameters is crucial for the validation of such a phenomenological model. 

In this paper, cornerstones of such a network model will be proposed. For an exploratory 
study a two-dimensional model suffices. It should give qualitative insight, but does not allow 
quantitative comparison with experiments, which will require a three-dimensional version 
and further refinements of the model in the future. 

The nanopowder is represented by a square grid, the nodes of which are partially filled 
with particles. On two opposite sides periodic boundary conditions are imposed. On the 
other two sides electrodes are attached to the powder. A certain electrical resistance is 
attributed to each bond connecting a particle to an electrode or to another particle. The top 
electrode can move towards the bottom electrode to simulate the effect of applied pressure. 

Possible mechanisms of densification during sintering are lattice diffusion, grain boundary 
diffusion, plastic deformation due to dislocation activity, sliding of particles, or viscous flow 
of completely or partially molten particles. The fact that a densification up to 97 % of the 
maximally achievable density is obtained within the short period of three minutes^ rules out 
the slow diffusion mechanisms for field assisted sintering. In a covalently bound material, 
the mobility of dislocations and grain boundaries can be neglected, too. Therefore, the 
most plausible model for the description of silicon nanoparticle sintering should be based on 
viscous flow. 

There are a couple of assumptions that may describe viscous flow during silicon nanopar- 
ticle densification: i) It has been demonstrated by differential scanning analysis (DSC) and 
in-situ transmission electron microscopy (TEM) that silicon nanopowders comparable in 
size and morphology partially melt and recrystallize upon heating at temperatures below 
1000 °C. It was observed, that such a partially molten powder could easily flow^. ii) A 
surface melting induced by the FAST process would result in a flow of particles™ 11 , iii) It 
is a plausible assumption that the silicon nanopowder is covered by a surface layer of native 
oxide. It is assumed that the amorphous oxide becomes viscous at elevated temperatures. 
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From the viewpoint of a model, it is of minor relevance, whether or not the surface layer 
is an amorphous oxide shell or native silicon. Both, molten silicon and softened silica would 
have a similar effect. Even if the material transport would happen by evaporation and 
redepositiorP^, the implementation into the theoretical model would be the same, so that 
this case is covered, too. 

For all lines of argumentation, the integration into a simulation model requires the same 
principal mechanism: Within the percolation path of the current, there will be a constellation 
of neighboring particles, where the ohmic heating will create the highest temperature. Above 
a certain threshold, this hottest nanoparticle will flow, due to complete melting, partial 
surface melting, or viscous flow of a surface layer. The matrix surrounding this flowing 
particle remains rigid, as the local temperature is lower there. 

From in situ TEM investigations 15 we assume that the flowing particle moves out of the 
current path. This also generally reflects the behavior of liquid droplets connected to solids 
(e.g. experienced when soldering). Within the real silicon nanocomposite, a recrystallization 
at adjacent particles with lower temperature will occur. If a native oxide shell is assumed 
to be the viscous phase, a segregation of the oxide accompanies this rearrangement. Within 
the simulation, we imagine that partially molten particles move into a neighboring void 
as demonstrated in figure |2j A stepwise rearrangement of the pores is the consequence. 
If a complete horizontal path is emptied, the upper electrode compacts the material by a 
movement downwards. 

At this stage of model development we do not implement heat diffusion explicitly. Instead, 
we assume that temperature differences within the sample, if not renewed by Joule heating, 
would vanish after a relaxation time r. This limits the buildup of temperature fluctuations 
in the sample. In the simulations we take r as discretization time. 

II. 2. Simulation Parameters 

In this section, the simulation model will be explained: 

We simulate the powder on a two dimensional square grid of size L x x L y = 200 x 250 
with periodic boundary condition in x-direction. The electrodes are the boundaries in y- 
direction. Each node of the grid can be empty (representing pore volume) or occupied by 
one particle. Hence the fixed grid spacing can be identified with the particle diameter. Size 
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differences between the particles are neglected. Particles sitting on nearest neighbor nodes 
(in x- or y-direction) are regarded as being in contact. 

The initial configuration, FigjlJ;) , is obtained by the following procedure: We start from 
a lattice gas, where half of the nodes are filled with particles at random positions, FigjT^i). 
The particles in contact to the upper electrode are shown in light blue. As long as there 
is no percolation path down to the lower electrode, the upper electrode can push the light 
blue cluster down by one grid spacing, without changing its structure. When this process is 
repeated, the cluster of particles in contact with the upper electrode grows (Figjl^)), until 
the first percolating path forms, which spans from the upper to the lower electrode (Figjlp)). 
It is assumed that the pressure applied to the upper electrode is so weak that a single 
percolating path is enough to resist further compaction. This is the initial configuration for 
the simulation of field assisted sintering. The initial inhomogeneity near the lower electrode 
is greatly reduced after a few simulation time steps as will be shown below. 

The simulation proceeds as follows: The current through the system is calculated by a fast 
algorithm described iiP^ for fixed voltage between the electrodes. The electrical resistance 
R attributed to the contacts between particles is taken as constant. The Joule heat RlfjT 
produced by the current from particle i to its neighbor particle j during one time step 
t is assumed to be delivered in equal parts to both particles, so that particle i receives the 
heat AQ i = 0.5 r. 

We assume that the the sample is kept at a constant average temperature T Sample and 
the heating compensates the energy loss to the environment at fixed T envmnt . 

1 N 

~jy ^ ] AQi = (AQ) oc r • (T samp i e — T envmnt ) (1) 

For T samp ie >> T cnvmnt we get the approximation 

(AQ) oc T sample (2) 

for the sample temperature. 

Particle % begins to melt, when the heat AQi exceeds the threshold 

Qthres. C ' (-^mclt -^sample) C ' -^sample ' I Tp ^ ) (3) 

\ -L sample / 

needed to raise the local temperature up to the melting point of the nanoparticles, T mc i t . 
The temperature dependence of heat capacity C close to melting has been neglected, and 



latent heat needs not be taken into account, since the onset of melting is the appropriate 
criterion for a particle displacement. 

With equation ([2l we get the melting criterion 



This form is computationally particularly convenient, because the ratio on the left hand 
side does not change, if the applied voltage is rescaled. This has the advantage that (AQ) 
does not need to be kept constant in eq. ([2]). Any change of (AQ) can be corrected by an 
appropriate rescaling of the applied voltage without any effect on the simulation. 

A particle with molten surface becomes mobile. It can be squeezed to any free neighboring 
place as shown in figure [2} This may interrupt the path of the electrical current. If more 
than one particle have molten they move in random order. Therefore it is possible that 
a disintegrated current path instantly reappears. After the movement, the particles are 
reintegrated into the rigid square grid. The grid has a lower mean temperature than the 
moving particles. Therefore, the reintegration physically corresponds to a recrystallizationP^l. 
We like to point out, that the model remains applicable, if instead of surface melting a 
temperature dependent viscosity of the surface layer leads to mobilization of hot particles 
and their subsequent arrest at a cooler neighboring place. 

If the electrical contact between the electrodes was lost because of particle movement, 
mechanical stability is lost as well. Then the cluster of particles in contact with the upper 
electrode is pushed down until percolation between the electrodes reappears. Hence the 
sample becomes more compact. 

Calculation of the currents, application of the melt criterion to move hot particles, and 
finally a compaction of the sample are repeated in every time step. The final configuration 
is reached, when the melt criterion is nowhere fulfilled any more or the sample became 
so dense that further compaction is impossible. Remarkably, the only model parameter is 



has so far not been sufficiently appreciated. Its effect on the final microstructure will be 
investigated in the following. 

As pointed out above, the initial configurations for our simulations have a zone near the 
lower electrode that is poorly connected to the upper one (figure [TJ^). This is a result of 
the method, by which the configurations were obtained, and will not be found in a real 




(4) 




The importance of this key parameter for field assisted sintering 
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powder. We characterize this inhomogeneity in figure |3j which shows the amount of free 
particles which have no contact to the upper electrode as function of the height (y-direction) 
for different stages of the simulation. As one can see, the inhomogeneity is strongly reduced 
after a short period of about 50 time steps, but it does not vanish completely. We assume 
that this has only a minor effect on the path formation, because after the first 50 time steps 
the cluster is connected to the bottom electrode over the width of the whole system, so that 
the path has enough freedom to form. 

We use a length unit of 20 nm for the grid to compare the simulation pictures with the 
experimental TEM images of the sample. This length corresponds to the average diameter 
of the nanoparticles that were used in the experiment. 

II. 3. Experimental Procedure 

Different batches of highly doped silicon nanoparticles with a nominal doping concentra- 
tion of 5 x 10~ 20 /cm 3 5 ■ 10 20 cm -3 and a mean starting diameter between 20 nm and 40 nm 
were used for field assisted sintering experiments, using a commercial system from FCT 
(Type HP D40, FCT Systeme Rauenstein, Germany). System and process have been ana- 
lyzed e.g. in Ref.^. Powders were precompacted and usually they have a density of around 
50 % of the value of crystalline silicon after this precompaction. Maximum temperatures 
were set between 860 °C and 1160°C at a heating rate of lOOKmin -1 , but it is known that 
the measurement technique by pyrometer tends to underestimate the actual temperature 
slightly^. Hold time was kept constant at 3 min. A pressure of 35 MPa was applied during 
the whole sintering process. Care was taken that the sample was electrically isolated from 
the graphite crucible. This was done by placing the sample within a graphite foil which was 
coated by boron nitride from the outside. The current was thus lead through the sample 
and not through the graphite crucible. The process was temperature controlled using a 
pyrometer focussed on the outer wall of the die wall surface. Current and voltage were set 
according to the time-temperature program, current being in the range of 0.5 kA to lkA, 
voltage in the range of a few V. 

Sintered nanocrystalline samples have about one third of the phosphorus atoms elec- 
trically activated, resulting in a charge carrier concentration around 1 x 10 20 cm -3 and a 
specific conductivity between 100 S cm -1 and 1000 Scm -1 . At room temperature, typical 
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values for the Seebeck coefficient are between 100 uVK -1 - 150 uVK -1 and for the thermal 
conductivity between lOW/m/K) - 20W/m/K). The high specific conductivity allows to 
conduct the experiment as described above. 

For an investigation of the microstructure, samples were cut into pieces and polished down 
to a thickness of around 40 pm. Final thinning to electron-transparency was done using a 
precision ion polishing system by Gatan inc (PIPS 691). Samples were characterized by 
transmission electron microscopy and scanning electron microscopy (SEM), the latter using 
a special adapter for TEM-samples. Microscopes used were a FEI Tecnai ST20 and Jeol 
JSM 7500F, respectively. The geometric density was obtained by the Archimedes principle 
by weighing the sample in ethanol and air. 

III. RESULTS AND DISCUSSION 

In this section we compare the results of the simulation with microstructural features 
found in the sintered experimental samples. Figure |4^i and [4j3 show an evaluation of the 
sintering protocol during densification of a nanocrystalline sample. A constant heating rate 
is followed by a hold time at 1060 °C. Two different phases can be distinguished: 

[I] The first phase is characterized by a constant density of the sample. A rearrangement 
of particles is highly unlikely as it would directly affect the porosity. Any changes of 
the electrical properties of the sample in this phase could be caused by thermal and / 
or electrical breakdown of the oxide layer^U^ and local welding at the micro-contacts 
due to Joule heating as found for larger particles^. 

[II] The second phase is characterized by a steep rise of density. Starting at the tempera- 
ture of 800 °C, the density increases from an initial value of 53 % up to 97 % caused 
by the sintering. This process slows down and the densification saturates eventually. 
The densification in this 2 nd phase is captured by the simulation assuming particle 
reorganization due to local melting. 

Around 750 °C to 800 °C, first melting events were found for very similar batches of silicon 
nanoparticles earlier^. Therefore, the temperature window fits well to a fluid phase assisted 
sintering model assuming partly or even completely molten silicon nanoparticles. 
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Figure |4]i shows the simulation results for the time evolution of the density for two 
different, constant sample temperatures. They are represented by the parameters m t h re s. 
1 2 and m t h re s. — 24 chosen such that the densification begins instantaneously. The density 
increases and saturates to a value that depends on the sample temperature, in agreement 
with phase II of the experiment. 

Figure [4^ ) shows the simulated final density as a function of m thres .. A low threshold 
represents a high sample temperature that is close to the threshold temperature for flowing. 
Accordingly higher thresholds are experimentally realized by lower sintering temperatures. 
Obviously, for high enough thresholds, no compaction takes place, in agreement with the 
low temperature behavior (phase I) of the experiment. Below a certain value of m t hres. 
the density starts to increase. Surprisingly, the simulation predicts that the compaction 
is maximal for a value of m t h res . ~ 12. For higher sample temperatures the final density 
decreases again. Experimentally it has so far only been seen that with increasing sample 
temperature also the compaction becomes stronger. Perhaps the high temperature behavior 
cannot be seen in the experiment, because the melting temperature is not uniform due to 
polydispersity. Another reason might be that coordinated movement of particles was not 
considered in the simulation, therefore a real material has more possibilities to close pores. 

The mechanism underlying the nonmonotonic temperature dependence of the final den- 
sity is the following: Let us consider cross sections of the simulated square grid parallel to 
the electrodes. Most heat is deposited in the cross section, which contains the least bonds, 
because this is the bottleneck of the current. For large m t hres., selectively bonds in this 
bottleneck melt until percolation between the electrodes is disrupted and the sample comp- 
tifies. This continues until the fluctuations in AQ t become too weak so that the threshold 
is no longer reached. Lowering m t h re s. makes it easier for compacted samples to still reach 
the threshold. However, this also widens the molten region, which no longer is restricted 
to the bottlenecks. This can be seen in Figure [5^,), where basically a whole current path 
melts. The corresponding particle displacement then amounts to a lateral diffusion of the 
percolating path rather than its disruption. Hence compaction is hindered. This explana- 
tion is supported by Figure [sJd) , where the accumulated heat is color coded to illustrate the 
diffusion of the path. 

While the development of the density seems to support our model, the change of the 
electrical properties of the sample during the ongoing FAST-process cannot be captured by 
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our simulation. A possible reason is that the moderate heating in phase I has not been 
included in the simulation. Within this phase, the percolation paths are not expected to 
disintegrate and fluctuate, as disintegration would be accompanied by a rearrangement of 
the powder network and thus a decrease of porosity which is not observed. Instead, existing 
paths 'burn' into the sample. Contact resistances within those paths may be optimized 
by breakthrough of the native oxide layer and by local welding. This optimization creates 
initial conditions for phase II which cannot be mimicked by constant contact resistances 
between the particles as we have assumed within the simulation. It can be concluded that 
the simulation so far qualitatively captures the compaction of the field assisted sintering of 
the highly doped silicon nanoparticles, but fails to model the resistivity decay and increase, 
because of the constant contact resistance. 

In the following paragraphs, we will analyze the microstructures obtained in field assisted 
samples. We used the contrast of backscattered electrons in SEM to make a fluctuation of 
the sample's density visible (figure [5^ and d), assuming that the atomic weight does not 
show a fluctuation^. Obviously, there are paths within the nanosilicon sample which have 
a higher density than the surrounding matrix, appearing at a micrometer length scale. In 
top view (figure [5] c) a honeycomb like structure is found. In the crosssection linear paths 
are found, which are very similar to the simulated heat pattern (figure [5] b). Indeed, in the 
simulation a weak density increase correlated to this heat pattern has been observed, Figure 
[5] c) . Such patterns are found regularly in nanosilicon samples produced in our laboratory. 
We conclude that these regions of higher density have their origin in the accumulated heat 
delivery (figure ^jp) produced by the electrical current which have 'burned' into the material. 
This accumulated heating over several time steps leads to a higher sintering activity of the 
hot regions and therefore to a higher local density. 

As a result, the simulation has proven to be in qualitative agreement with experiment, 
as far as the developing microstructure is concerned: Accumulated heat fluctuations are 
correctly given by the model and found as density fluctuations in the experiment. 

In the following, the parameters of the simulation are varied and compared with ex- 
perimental investigations. Therefore, the final microstructure of the simulation and the 
experiment are compared in figure |6j Two threshold-values have been chosen, m t h re s. = 12 
and m t hres. = 24, to represent a medium sintering temperature (close to optimal compaction) 
and a low one, respectively. The medium sintering temperature resulted in a relatively ho- 
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mogeneous sample, both in simulation (figure |6pi) and experiment (figure [6p and c). Figure 
[6J3 shows a nanosilicon which was sintered at 1050 °C in transmission. The nanocrystalline 
structure is preserved despite of the density fluctuations due to the current paths as discussed 
above. 

For the low sintering temperature, a different microstructure is expected according to the 
simulation. Regions which have been molten and regions which were not heated sufficiently 
are found adjacent to each other and seem to cluster as seen in figure [6ji. This is found 
by a SEM analysis of the accordingly sintered nanosilicon. On a scale of micrometer, very 
well-sintered as well as bad-sintered regions have formed (figure [6]f ) . 

The reason for this may be attributed to the temperature gradients created within the 
sample. If we assume that the viscous flow of the particle physically always starts at the same 
threshold temperature (e.g. the melting temperature of the nanoparticles), a low sintering 
temperature (or high mthres.) means that the temperature variation within the sample is 
higher because the mean temperature is lower. This will necessarily create inhomogeneities. 
Note, that these inhomogeneities cannot be reduced by longer sintering times, nor will the 
density be increased considerably. The well conducting paths have burned into the sample 
after a short time. The bad- sintered regions will not carry the current to the same amount 
as the well-sintered regions. Thus, longer sintering times may even result in more extreme 
inhomogeneities. It can be concluded: If the mechanism of densification is viscous flow, 
homogeneous samples can only be obtained by choosing the sintering temperature close to 
the threshold temperature of the viscous flowing. 

Figure [7] shows a detailed analysis of the temperature distribution (number of nodes n 
at heat AQi) during the simulated sintering for m t hres. = 12 and m t hres. = 24. It can be 
seen that the temperature distribution for m t hres. = 12 remains very close to the initial 
configuration. Opposing, the distribution changes during the process for m thres = 24. More 
nodes with lower AQi are found and less particles become mobile during each time step. 
What does this mean for field assisted sintering? There are only few hot spots created during 
the sintering if the mean temperature is considerably lower than the threshold for flowing 
(m t hres. = 24). At these hot spots, the microstructure will be reorganized only very locally, 
leading to well-sintered areas which in turn will carry the current afterwards. 

For a homogeneously sintered sample, it is important that current paths not only form, 
but disintegrate again. If the mean temperature during the experiment is close to the thresh- 
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old for flowing (m t hres. = 12), many particles become mobile in each time step. Initially con- 
ducting current paths will thus be destroyed by the movement of the flowing particles. That 
is good for the formation of a homogeneous microstructure. Of course, if a nanocrystalline 
material is intended, the temperature profile has to be well balanced: Too many flowing 
particles will lead to an unwanted coarsening. 

IV. CONCLUSION 

Concluding, we used a network model to simulate the field assisted sintering of highly 
doped silicon nanopowder and compared the simulated results with experimental data. The 
evolution of density during the sintering process qualitatively matched. We found a good 
agreement between experiment and simulation with respect to the local temperature distri- 
bution which is displayed as density fluctuation in the microstructure of densely sintered 
nanocrystalline silicon. So we captured some of the main aspects of the FAST process with 
our simulation approach. As a consequence, some information can be extracted as a guide 
for sintering dense, homogeneous bulk applying field assisted liquid phase sintering: After a 
very fast sintering phase in which the densification is observed, the density saturates. Longer 
hold times, especially at temperatures considerably below the threshold for flowing, will not 
lead to a further densification but rather increase the microstructural inhomogeneities within 
the sample. For the developing microstructure, short sintering times at elevated tempera- 
tures — i.e. a temperature close to the threshold for flowing — result in best homogeneity 
combined with high density. 
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FIGURES 




FIG. 1. Two steps of the initial compaction until the powder percolates through the system. Empty 
nodes are colored black, the electrodes are green. The light blue nodes belong to the moving cluster 
with contact to the upper electrode. The dark blue nodes are occupied. The resulting percolating 
network is shown in picture c) in cyan. 




FIG. 2. Possible movements of a flowing particle. The color coding gives a rough example of the 
temperature of the particle at the given time step. White squares represent molten nodes/particles, 
yellow color indicates nodes which have been heated up and black color represents cold nodes. 
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FIG. 3. Free particles, e.g not connected to the spanning percolation cluster, versus the distance 
from the bottom electrode. The different curves are after different number of time steps of the 
sinter model algorithm for m t hres. = 12. Curve 1 shows the initial situation. The inhomogeneity 
caused by the initial compaction is obvious. Curve 2 shows the situation after 50 time steps, curve 
3 after 400 time steps. The initial increase of free particles towards the bottom electrode is strongly 
reduced already after 50 time steps. 
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Experiment 



Simulation 




FIG. 4. (a) Temperature T and Joule heat Q delivered per unit time, vs. time. T is measured 
by a pyrometer focused on the upper graphite die. Q is calculated from the applied current 
and voltage. The shown process is characterized by a constant heating rate of lOOKmin^ 1 and 
a 3min hold time at 1060 °C. (b) Density of the sample. The density is calculated from the 
movement of the upper graphite die and related to the end porosity derived from the geometrical 
density by Archimedes principal. Different phases of the process can be distinguished. See text 
for explanation, (c) Theoretical final densities received from simulations for different m t hres.- (d) 
Densities development received from two different simulations with m-thres. = 12, and m-thres. = 24 
for the densification (Phase II) of b). 
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FIG. 5. (a) Heat AQi generated during the 100 time step of a simulated sintering process for 
^thres. = 5. Particles plotted in green are above the threshold, have molten and will be moved as 
indicated in the text, (b) Accumulated heat at each particle after 400 time steps for for mthres. = 5. 
(c) Density in a vertical slide through the system shown in b). To reduce the noise, the slide has 
a width of 1.6pm, which is in the same order as the path observable in b). (d) SEM image of 
a nanosilicon sample compacted by field assisted sintering (1060 °C, 3min). The same sample is 
shown as cross-section in e). Backscattered electrons are used for imaging to indicate fluctuations 
in density. Denser regions appear brighter. 
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FIG. 6. (a) Powder compacted with mthres. = 12 which corresponds to a high sintering tem- 
perature. The red region consists of molten particles, while the grey dots are particles that are 
still in their original condition, (b) Scanning TEM (STEM) image of a nanosilicon with sintering 
temperature of 1060 °C, density 97.2 %. In STEM mode it can be seen, that the nanosilicon sample 
still has a high porosity, (c) Besides, SEM investigation shows that the samples has a good overall 
homogeneity — inclusions seen in this image appear to be areas enriched with oxygen, (d) Powder 
compacted with mthres. = 24 which corresponds to a low sintering temperature. The color-coding 
is the same as in picture (a). The sample has a much higher amount of not-sintered powder (grey). 
Molten regions and not-sintered powder seem to cluster, (e) STEM image of a nanosilicon with 
sintering temperature of 860 °C (density 95.0 %), demonstrating porosity, (f) A SEM investiga- 
tion showed that the sample is very inhomogeneous on the micrometer scale. Large, well sintered 
regions have formed and are found adjacent to regions of not-sintered powder. 
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FIG. 7. Number of nodes with a normalized heat AQj. Curve 1 represents the initial distribution. 
The initial distribution shows no dependency of m t hres. and is only displayed for m t hres. = 12 . 
Curve 2 shows the heat distribution after 400 pulses for m t hres. = 12, curve 3 the same situation 
for m t hres. = 24. While there is not much change in the heat distribution of the sample in the case 
of a low threshold (mthres. = 12), it is found that the overall heat distribution changes drastically 
for higher thresholds (mthres. = 24). Only few particles are mobile (in their liquid phase), while a 
much higher number of nodes have a relatively low AQj, creating an inhomogeneous temperature 
field within the sample. 
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